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Abstract 

A FORM based package (ON-SHELL2) for the calculation of two loop self-energy 
diagrams with one nonzero mass in internal lines and the external momentum on 
the same mass shell is elaborated. The algorithm, based on recurrence relations 
obtained from the integration-by-parts method, allows us to reduce diagrams with 
arbitrary indices (powers of scalar propagators) to a set of master integrals. The 
SHELL2 package is used for the calculation of special types of diagrams. Analytical 
results for master integrals are collected. 
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1 Introduction 



High experimental accuracy achieved in the last years allows to test the Standard Model on 
the level of quantum corrections. Therefore to improve the accuracy of predictions within 
this model is of urgent need. The calculation of mass- dependent radiative corrections is 
complicated but can be performed to a large extent by using computer algebra. In recent 
years many algorithms have been developed and huge program packages were elaborated 
for this purpose (for a review of existing packages see Ref. [1]). Due to different mass 
scales of the particles in the Standard Model the method of asymptotic expansion [2] 
- the expansion of Feynman diagrams w.r.t. the ratio of different scale parameters - is 
becoming more and more popular. The calculation of radiative corrections to low energy 
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processes e.g., in particular those with light external fermions, in many cases reduces 
to the calculation of self-energy diagrams with external momentum at different scales. 
This is one of the reasons why the evaluation of two-loop self-energy diagrams is worth 
special attention. From the point of view of approximation methods two-loop self-energy 
diagrams can be divided into several classes: 

• Only one non-zero mass enters internal lines and the external momentum is on the same 
mass shell. The calculation of diagrams of this type occurring in QED and QCD has 
been implemented as the package SHELL2 [5] . 

• There are heavy particles in internal lines and the external momentum is on the mass- 
shell of a light particle [6-8]. Diagrams of this type occurring in the Standard Model 
have been collected in the package TLAMM [9]. 

• All internal particles are light or massless and the external momentum is on a heavy 
mass shell (see, for example, Ref. [10]). 

• Several different heavy masses occur in internal lines and the external momentum is on 
the mass shell of a heavy particle [11,12]. 

• Diagrams close to thresholds (see Ref. [13] and references therein). 

A general recipe of reduction of arbitrary two-loop self-energy diagrams to a set of master 
integrals has been suggested by Tarasov [14]. The algorithm was implemented in FORM 
and then in MATHEMATICA [15] . However, this method suffers from the drawback 
that in processing the reduction of scalar master integrals with shifted dimension to the 
generic dimension of space-time, powers of - may arise which require the expansion of 
master integrals as series in e, which is a difficult task. This problem is avoided in our 
approach. 

We present a FORM [16] based package that allows to calculate arbitrary two-loop self 
energy diagrams with one non-zero mass and the external momentum on the (nonvanish- 
ing) mass shell. Our algorithm concerning the V-type diagrams is very similar to the one 
described in Ref. [17]. 

The paper is organized as follows. In Sect. 2 the full set of recurrence relations is presented 
which allows to express the initial diagrams with scalar product in the numerator in terms 
of diagrams with positive (> 1) indices. Sect. 3 is devoted to the description of how to 
use the package. In appendix A (appendix B) we give all needed recurrence relations to 
reduce the scalar F-prototypes (V-prototypes) with arbitrary positive indices to a set of 
master- integrals. In appendix C the analytical results for all integrals shown in Fig.l with 
indices 1 are collected. Even though not all of these are considered as master integrals, 
they can be used for comparison. We are working in Euclidean space-time with dimension 
N = A — 2e. 



4 One of the first calculations of this type for the Standard Model and QED was performed in 
Refs.[3,4]. 
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2 The recurrence relations. 
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Fig. 1. The F, V and J topologies. Bold and thin lines correspond to the mass and massless 
propagators, respectively. 

The full set of two-loop self energy diagrams with one mass and external momentum 
on the same mass shell is given in Fig.l. We distinguish three basic topologies which in 
accordance with notations in Ref. [14] we call F, V and J prototypes with five, four and 
three lines, respectively. Our notation is given in Fig. 2. The diagrams implemented in 
the package SHELL2 (F01101, F00110, V1110, VOOll, V1000 in our notation) and those 
considered in detail in Refs.[18]-[20] (FOOOOO, VOOOO, JOOl, JOOO) are not discussed here 

n 

5 The procedures for the calculation of all diagrams of the topologies shown in Fig.l are imple- 
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Fig. 2. Notations used in present paper. 



The general prototype involves arbitrary integer powers of the scalar denominators cl = 
powers 2l are called indices of the lines. The mass-shell condition 
for the external momentum now is p 2 = —m 2 . Any scalar products of the momenta in 
the numerator arising from projection or expansion are reduced to powers of the scalar 
propagators (in case of V and J topologies the corresponding lines are added). Thus, 
the indices may sometimes become negative. Recurrence relations are derived via the 
integration- by-parts method [19] and applied to the massive case as in Ref. [22]. They 
allow to reduce all lines with negative indices to zero and the positive indices to one or 
zero. Further we use the shorthand notation {123} of Ref. [23] to denote the relation for 
the triangle formed of lines 1, 2, and 3: 



d N k / , 2m? m? + ml — m? 9 + Ci 2 — Ci 



c 2 



m\ + ml — m| 3 + c 13 - c% 
c 3 



{123} 



where a double index like {12} refers to a line that starts at the point where lines 1 and 
2 meet (see Fig. 3). For an external line on the mass shell, the value of cl is equal to zero. 



Ci, m 



c 12 , m 12 




Fig. 3. "Triangle" rule. 



mented in our package. 

6 Their explicit expressions are c\ = k\ + mf , c 2 = h\ + m|, C3 = (k\ — p) 2 + m 2 , C4 
(k 2 - p) 2 + ml, c 5 = (h - k 2 ) 2 + m| 
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2. 1 F -topology 



To exclude the numerator (for example, ji < 0, c^ 1 in the numerator) of F-type diagrams 
we use the following set of recurrence relations: 



{245} ^ci = -^2m\ - J -± (m\ + m\ - m 2 - c 2 ) 
c 5 c 2 c 4 v - 1 



( m 5 + m 2 - m l ~ C 2) - N + 2j 2 + ji + j 5 , 



(2) 32 + 1 



{524} — ci = 2m 5 [m 5 + m 4 — m 3 + c 3 — c 5 

c 2 c 5 c 4 v 



(3) js ^ 1 



(777,5 + ml-ml- c 5 ) - N + 2j 5 + j 4 + j 2 , 
c 2 v y 



{245} + {135} — ci = — (ml + m\ - m 2 ) 

c 3 c 3 

i 3 4 / 2 1 2 2\ 

H — (^m 2 + m 4 — m J 

_!± C2 + Zl 2 m? + ^2m^ + ^2m* 
C 4 Ci c 2 c 5 

+2N-2j 1 -2j 2 -j 3 -j A -2j 5 , 



(4) j 2 = js = k = 1 



{315} iV - 2j 3 - Ji - 3n = ~- (ml + m\ - m 2 - c 3 ) 

^5/2, 2 2, \ ^3 n 2 

( m 3 + m 5 — m 4 + c 4 — c 3 2m 3 , 

C5 v ' 7 c 3 

where both sides of these relations are understood to be multiplied by 

/ ifi^huh ■ The relations for (32,33,34) < are obtained from symmetry properties of 

c l c 2 c 3 c 4 c 5 

the integral under consideration: 



(ji, m i) <-»■ (Js, m 3 ), (j 2 , m 2 ) <-> (j 4 , m 4 ), 

and 

(ji, mi) <-> (j 2 , m 2 ), (j 3 , m 3 ) <-> (j 4 , m 4 ). 
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C5 from the numerator can be eliminated by a general projection-operator method [19]. 
Using the decomposition 

kik 2 = A(k u k 2) p) H - 9 , 

where A(ki, k 2 ,p) = k% (d^ — k 2 , and the property that odd powers of A(ki, k 2 ,p) 
drop out after integration and for even powers we have 



J d N h d N k 2 M^p] f 2 [k 2 ,p] A 2n (k u k 2 ,p) 
r (n + \) T[\(N - 1)] 

;§)r[n + ±(jv-i)] i=1 . 



r(n + i)r[i(iv-i)] i r , ,*»,,, \ 



it is possible to reduce the initial integral to a product of one-loop integrals. 

Using the above relations, F-type integrals with arbitrary indices are reduced to F-type 
integrals with only positive indices or V-type integrals with arbitrary indices. For the 
former case a proper arrangement of recurrence relations in general reduces the sum of 
all indices by 1. These relations are given in appendix A. Only eight diagrams (Flllll, 
F00111, F10101, F10110, F01100, F00101, F10100, F00001) with all indices equal 
to 1 form the basis for F-type diagrams. 



2.2 V-topology 



Consider now the V-type diagrams. The recurrence relations we use are {425}, {423} and 
the following set: 



{530} N - 2j 5 - j 3 + -1m\ + 3 - (m 2 3 + m\ - m\ + c 4 - c 5 ) = 0, 



J -±2m 2 3 + * 
C3 c 5 

h ( m l - m\ + m 2 + c 2 — c 4 ) (ml - m\ - ml + c 3 + c 4 - c 5 ) 
c 5 2c 6 

+ il 2m 2 - ( * + ?± + J l) (ml _ ml + m 2 + C2 _ \ = 

c 2 \c 2 c 4 c 5 / v " y 



{350} N - 2j 3 - j 5 + 3 -^2m 2 3 + ^ (m* + m\ - m\ + c 4 - c 3 ) = 0, 

C3 c 5 v 1 7 

„• 1^2 _2 1 

{^} 



where c@ = c 4 — m\ (see Ref. [17]). If 7^ 0, the expression can be simplified later 
by partial fraction decomposition. For all cases (j 2 , j 3 , j 5 ) < 0, except j 4 < 0, the initial 
diagram can be reduced to two-loop tadpole-like integrals by means of [6]: 
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d N k 2 (k 2 2 ) j f(k 2 ,k 1 ) 



For j 4 < we write c 4 = c 4 + m\ and redefine c 4 — c 4 , which allows to consider only the 
massless case. Then the following recurrence relations are needed: 



{350} ^c 4 = 3 -lU-ml- O - 2-^3 ~N + 2j 3 + 35, 

(2) J3 7^ 1 



{530} ^c 4 = ^ (c 5 - ml - ml) - 2^m\ -N + 2j 5 + j 3 , 
C3 c 3 v 1 c 5 



(3) J2 7^ 1 



{425} + {350} — c 4 = — 2m\ + — 2m l + — ^\ 

c 2 c 3 c 4 c 5 



+ ^ ( m 2 _ m A + 2iv - j 2 - 2j 3 - 2j 4 - 2j 5 . 

(4) j 2 = j 3 = h = 1 



25 + 2{425} + {350} 3iV - 4j 2 - 2j 3 - 2j 4 - 2j 5 = -2ml— ~ 4m 2 2 — 

c 3 c 2 

, 35 C2 / |2 2\ , J5 C 3 / 2 2 \ 

H C4-C3 + 1T1, -m. H m 9 — m 

c 5 c 4 V 3 5 J c 5 c 4 V V 

+ U5 + 2j 4 ) h — (m -m 2 - 2m 5 

c 4 c 5 v 

| j 5 (m| - roj) (m 2 - m|) 



c 5 c 4 

The result of application of the above recurrence relations are V-type diagrams with 
only positive indices or J-type integrals with arbitrary indices. The full set of recurrence 
relations for the former case is given in appendix B. The complete set of basic integrals 
is just given by Villi and V1001 with indices equal to I. 



2.3 J -topology 



The integrals of this type are discussed in detail in Refs. [13,14]. We only mention here, 
that to reduce the numerator the following recurrence relation, suggested by Tarasov [24], 
is needed: 
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(N + U! + U2 -2)v{v 1 ,v 2 ) =P 2 {{vi - l)k 2 1 l-+iy 2 (k 1 k 2 )2 }1 v{v u v 2 ), 



where v(vi,v 2 ) = j d N kid N k 2 f(k 1 ,k 2 )(k 1 p) Ul (k 2 p) l/2 and f(k 1: k 2 ) is an arbitrary scalar 
function; v\, v 2 > and l ± t>(z/i, u 2 ) = t>(z/i ± 1, z/ 2 ), etc. 

The master integrals are the following: one prototype Jill with all indices equal to 1, 
and two integrals of JOll-type: with indices 111 and 112, respectively. 

2.4 Master-integrals 

To obtain the finite part of two-loop physical results one needs to know the finite part of 
F-type integrals, V- and J-type integrals up to order e, and one-loop integrals up to order 
e 2 . A detailed discussion of the calculation of master-integrals is given in [21]. Here me 
mention only, that the calculation of the e (e 2 ) parts has been performed by the differential 
equation method [22]. The results are collected in Appendix C. 



3 Use of the package 

The package consists of a set of procedures for the calculation of all two-loop integrals, 
presented in Fig.l (flllll.prc, on3.prc, on2.prc, etc), two-loop tadpoles (villi. pre, 
vlOll.prc, vlOOl.prc) and one-loop integrals (vll.prc, onl.prc, onsll.prc), where "1"("0") 
in the name of the procedure stands for massive (massless) lines, respectively. on3 and 
on2 are two-loop integrals from the SHELL2 package, vll is the one-loop massive bubble, 
onsll and onl denote the one-loop self-energy on-shell integrals with two and one massive 
lines, respectively. The integration momenta in the package are denoted by Kl and K2 
for two-loop integrals and by Kl for one-loop integrals, P is the external momentum. All 
scalar products in the initial diagram must be rewritten in terms of propagators: 

ci — C3 + m\ — m\ — m 2 
2 ' 
c 2 — C4 + 777.4 — rn\ — rn 2 
2 ' 
c i + c 2 _ c 5 + m\ — m\ — m\ 
2 ' 

ci - ml, 
c 2 - m\. 

To specify the type of two- (one-) loop diagrams, the products of scalar propagators must 
be substituted by the proper functions of F-, V-, J- and ON-type with arguments denoting 



hp = 

k 2 p = 

u. u. — 
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the indices and a symbol for the mass shell. 



To work with fractions of N-dimensional numbers, two functions, SS and NN (originating 
from the package "LEO" [23] ) are used: 



NN(a,j) = (N + ay , j>0, 

and 



After application of each recurrence relation the procedure "ration" for the simplification 
of products of SS's and NN's must be called. The procedure "finitem" substitutes the 
values of master integrals and performs the expansion of the functions NN and SS in e. 

The integration procedure starts with F-type integrals. We apply the recurrence relations 
given explicitly in appendix A. After applying them several times, the integrand is reduced 
to the master integrals or to new, more simple, integrals like V-type, e.g. Then we apply 
several times the recurrence relations of Appendix B for the V-types. One needs to call 
all procedures step by step to reduce the initial diagram to the set of master integrals. 
The recommended sequence for calling the procedures is the following one: 

#call f 11111-C'TIMES'} 
#call fOllll-C'TIMES'} 
#call f 11110-C' TIMES'} 
#call f00111{' TIMES'} 
#call f 10101-C' TIMES'} 
#call f 10110-C' TIMES'} 
#call f01100{' TIMES'} 
#call f00101{' TIMES'} 
#call f 10100-C' TIMES'} 
#call fOOlOO-C' TIMES'} 
#call f00001{' TIMES'} 
#call f00000{' TIMES'} 
#call vllll-C' TIMES'} 
#call vOlll-C' TIMES'} 
#call vlOll-C' TIMES'} 
#call vl01(K' TIMES'} 
#call vOllO-C' TIMES'} 
#call vlOOl-C' TIMES'} 
#call vOOlO-C' TIMES'} 
#call vOOOl-C' TIMES'} 
#call vOOOO-C' TIMES'} 
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#call j 01 1{ 'TIMES'} 
#call on3{' TIMES'} 
#call on2{' TIMES'} 
#call vllll{' TIMES'} 
#call vl011{' TIMES'} 
#call vlOOl-C' TIMES'} 
#call onl{2*' TIMES'} 
#call onsll{2*' TIMES'} 
#call vlOO{' TIMES'} 
#call vll{2*' TIMES'} 

All programs of the package are realized with the help of FORM procedure facilities. To 
perform the integration, one needs to call the procedures with the name of the correspond- 
ing prototypes and one argument which determines how often the recurrence relations are 
to be called. This number depends on the complexity of the calculated diagram. In most 
cases it is equal to the sum of indices of the integrand. 




Example 1 
Fig. 4. 

Let us consider as an example, the physical diagram shown in Fig. 4 in more detail. The 
FORM input is created automatically by DIANA [25] . For two- loop self-energy diagrams 
DIANA generates all necessary information, e.g. identifying symbols for the particles 
of the diagram and their masses , distribution of integration momenta, Feynman rules 
(linear/nonlinear gauges), number of fermion loops, symmetry factors, etc. To calculate 
the transverse part it is sufficient to make the following substitutions: 

multiply, (d_(mu,nu)-p(mu)*p(nu)/p.p)*SS(-l, 1) ; 
. sort 

id kl.p = (kl.kl-c3)/2; 

id k2.p = (k2.k2-c4+mmZ-mmW)/2; 

id kl.k2 = -(c5 - kl.kl - k2.k2 - mmW)/2; 

id kl .kl = cl - mmZ; 

id k2.k2 = c2 - mmW; 

id p.p~j? = (-mmW)~j; 

. sort 

id mmZ~j? = mmW" j ; 
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id l/cl~jl?/c2~j2?/c3~j3?/c4~j4?/c5~j5? = FllllKjl , j2, j3, j4, j5,mmW) , 



where we have explicitly set m 2 ^ = m| Calling the above routines (not all of them are 
needed in this special case) yields the result 



Examplel = m^Flllll(l, 1,1,1, 1, m 2 w ) (33 - ' ' 



N- 1 



+V1111(1, 1, 1, 1, m 2 w ) -42 + ^—^ + 



2 27 



N - 1 4(iV- l) 2 y 



26 17 



m 2 w V 3 3(iV-l) 4(JV-l) 2 y 

| vmijui^ / 5 9 



+ 



m 2 w \ 2{N-l) A(N-l) 2 / 

2/5 21\ 



ONSll(l,l,m^ 



{N-l) 2 



+ ONSll(l, 1, mf„) (-3^ - ^ + 3^ - 
( 64 24 32 8 32 \ 



k 3(N-4) (iV-4) 2 (iV-2) (iV-2) 2 3(iV-l)/' 

4 2 

where the overall factor is ^g"^ ■ After calling "finitem" we have: 



1417 1 /10375 665 vr \ 21187 21 . . 88 >/N 

Examplel = -—r - - — — — -p — — + —( 2 - — C(3) 

F 24e 2 £ V 48 12 V3 J 32 8 s w 3 s w 

4007 7T 665 7T , 16449 „ „ nrt tt „ 

+ f= f= n 3 + S 2 + 132^S 2 

18 V3 12 y/3 16 v/3 

59.0 115.6 „ rt „ 
69.6. 



e 2 



So far we have considered only the case of merely one non-zero mass. Of course it is 
obvious that as further application of our package we use it for the expansion of diagrams 
in terms of mass differences. In general a 'standard' expansion of the scalar propagators 
in terms of the mass difference, i.e. in the above case in terms of m 2 ^ — m|, yields as 
expansion coefficients again integrals, which can be handled by our package. In order to 
demonstrate this possibility, we extend the calculation of the diagram in Fig. 4 up to the 
second order in A = 1 — m^/ml = sin 2 9w with the result 



1 / 1417 667 . 95 A ... 
Examplel = — ( — + —A - —A 2 
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1 /10375 24197 A 17957 2 \ 1 vr /665 137 A 635 2 \ 

7T /4007 21881 A 45067 a2 \ tt , /665 137 . 635 A 2 N 

A + A 2 = In 3 A A 2 

V3 V 18 72 432 J yft \ 12 3 72 , 

n ( 16449 7461 . 2195 A 2 \ ?r / 351 A 699 . 2 \ 

/21187 216049 A 403637 A 2 \ ^, n j21 1477 A 23179 . 9 \ 
-A + — — A 2 +C(2) + 



32 216 864 / ' v ' \ 8 36 432 

59 
2 



<(3) ( ^ - ?A + ^A 2 ^) » 1 (-59.0 + 83.4A - 23.8A 2 ) 



+- (-115.6 + 253.2A - 140.7A 2 ) - 69.6 + 211.6A - 118.4A 2 . 



With the numerical value A ~ 0.23, we see that the convergence in A of the above series 
in the three contributions is quite good and it can be expected that it will even improve 
due to 'gauge cancellations' if a complete gauge invariant subset of diagrams is taken into 
account. 



4 Conclusion 



The presented package has been developed for the calculation of two-loop self energy 
diagrams with only one non-zero mass. All mass combinations of the Standard Model 
with heavy masses are included. In this sense it is an extension of the existing package 
SHELL2, which takes into account only diagrams occurring in QED and QCD. We have 
also shown that this package can be used for the case of different masses in the SM 
by expanding in terms of mass differences. Thus we have provided a program for the 
evaluation of at least a large class of on-shell two-loop self energy diagrams in the SM. 
The time of calculation of one diagram, depending of course on the order of expansion in 
the mass differences, is not very large in general. 



Acknowledgments We are grateful to A. Davydychev, A. V. Kotikov, O. V. Tarasov, 
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Appendix 



A The set of recurrence relations for F-type integrals 



The full set of recurrence relations valid for arbitrary masses and momenta is given in 
Ref. [26]. Here we present some rearrangement of these relations which allows to reduce 
the positive indices of the on-shell diagrams with one mass to or 1. 

(1) Flllll 

4m 2 — = — (c 2 - c 3 - c 5 ) + — (3ci + c 4 - c 5 ) 

Ci Ci c 3 

+ - (3d - 3c 2 + c 4 - c 3 ) - N + 4ji, 

4m 2 — = — (3c 5 - 3c 2 - c 3 ) + — (3c 5 - 3c 4 - ci) 
c 5 Ci c 3 

+ - (c 4 - c 3 + c 2 - Cl ) - AT + 4j 5 . 

c 5 



(2) F01111 



iV - 2ji - j 3 -j 5 = — (ci - c 2 ) + — ci. 

c 5 c 3 



(3) F11110 



iV - 2j 5 - ji - j 3 = — (c 5 - c 2 ) + — (c 5 - c 4 ) 

Cl c 3 



(4) F00111 



2m 2 — = — (2c 5 - 2c 2 - c 3 ) + — (2c 5 - 2c 4 - 3ci) 

Cl Ci c 3 

+ - (3c 2 - 3ci + c 4 - c 3 ) + 2iV - 3j 3 - 5ji, 
c 5 

o 2 J3 Jl J3 . J5 i . \ . . 

2m — = — c 3 ci H (c 2 - ci + c 3 - c 4 ) + j 3 - j u 

C3 Ci C 3 C 5 

2m 2^i = J3 Ci + J± C2 -2N + 2 3l + 2j 2 + 2j 5 + j 3 + j 4 . 

c 5 c 3 c 4 



(5) F10101 



2Jl J3 , v . Jl , \ 

m — = — (c 5 - c 4 ) H (c 5 - c 2 - c 3 J 

Cl c 3 Ci 
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+ — (C 4 - C 3 ) + j 5 - J 3 , 

m 2J2 = J5 ^ _ ^ _ J_2^ + ^ _ g _ ^ _ . 

C 2 C 5 C 2 

2m 2 — = — (c 5 - ci) + — (c 5 - c 3 ) - N + 2j 5 + j 2 + j 4 . 

c 5 c 2 c 4 



(6) F10110 



2m 2 — = — c 4 (c 3 - c 4 ) - iV + 2j 4 + j 2 + j 5 , 

C 4 C 2 C5 

"i 2 — = — (ci - c 5 ) + — (c 3 - c 5 ) + iV - 2j 5 - j 2 - j 4 , 

c 2 c 2 c 4 

ra 2 — = — (ci - c 2 ) - — c 2 + A^ — 2j 2 - j 4 - J 5 , 

c 5 c 5 c 4 

m 2^i = Jl _ ^ + J3 ^ _ ^ _ N + ^ + ^ + ^ 
Ci Ci c 3 

2m 2 — = — (c 4 - c 5 ) + — (c 3 - c 4 ) + — (c 2 + c 3 - c 5 ) + j 3 - j 5 

c 3 C 3 C 5 Ci 



(7) FOllOO 



2m 2 ^ = J - (2c 2 - 2c 5 + c 3 ) + - (2c 4 - 2c 5 + a) 

Ci Ci c 3 

+ — (ci - c 2 + c 3 - c 4 ) + ji + j 3 - 2j 5 , 

O 2^2 j 2 J4 J5 / . x 

2m — = — c 4 H c 2 H (c 2 - Ci + c 4 - c 3 ) 

c 2 c 2 c 4 c 5 

-2AT + 2j 5 + 3j 2 + 3j 4 , 

m 2 — = — (c 2 - ci) - — ci + N- 2ji - j 3 - j 5 - 

C5 c 5 c 3 



(8) F00101 



m 2 — = — (ci - c 2 ) + — ci - iV + 2ji + j 3 + j 5 , 

C5 c 5 c 3 

m 2 — = — (c 3 - c 4 ) - — c 4 + A^ — 2j 4 - j 2 - J5, 

C2 C 5 C 2 
2^4 J5 / x 3a .33 . n ■ O ' 1 ' 

m — = 2— (ci - c 2 ) c 2 H ci + 2ji - 2j 2 + j 3 - j 4 , 

c 4 c 5 c 4 c 3 

2m 2 — = — c 3 - 2— ci + — (2c 2 - 2ci + c 3 - c 4 ) + - 3ji - j 5 , 

c 3 Ci c 3 c 5 

m 2 — = — (c 5 - c 2 - c 3 ) + — (c 5 - c 4 ) + — (c 4 - c 3 ) + j 5 - j 3 . 

Ci Ci c 3 c 5 



(9) F10100 
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2m 2 — = — (c 5 - c 2 - c 3 ) + — (ci - c 4 + c 5 ) 
Ci Ci c 3 

+ - (ci -c 2 -c 3 + c 4 ) - N + 2ji + 2j 5 , 
^5 

2m 2 — = — (ci - c 4 - c 5 ) + — (c 2 + c 3 - c 5 ) + — (c 2 - ci + c 3 - c 4 ) + iV - 2j 4 

c 2 c 2 c 4 c 5 

2m 2 ^ = ^ (3c 2 + c 3 - 3c 5 ) + - (c 4 + 3c 4 - 3c 5 ) 

C 5 Ci c 3 

+ - (ci - c 2 + c 3 - c 4 ) + iV - 4j 5 . 

c 5 



(10) FOOIOO 



N - 2ji - j 3 - J 5 = — (ci - c 2 ) + — Ci. 

C5 c 3 



11) FOOOOl 



4m 2 — = — (c 5 - c 2 - 3c 3 ) + — (ci - c 4 + c 5 ) 

Ci Ci c 3 

+ - (ci - c 2 - 3c 3 + 3c 4 ) + iV - 4j 3 , 
^5 

4m 2 — = — (c 5 + c 3 - c 2 ) + — (c 5 + ci - c 4 ) 

c 5 Ci c 3 

+ - (ci - c 2 + c 3 - c 4 ) - 3N + 4ji + 4j 3 + 4j 5 . 

c 5 



B The set of recurrence relations for V-type integrals 



The full set of recurrence relations for the V-type integrals is also given in Ref.[26]. Here 
we give our rearrangement of these relations which again allows to reduce positive indices 
of the on-shell diagrams with one mass to or 1. 

(1) Villi 

3m 2 — = — (c 4 - c 5 ) + 2— (c 3 - c 4 ) - N + 3j 3 , 



c 3 c 3 c 5 



2J2 ■ C2 C 2 C 3 \ ? 4 75 J 2 . 

3m — = j 5 + 2— c 2 H c 2 c 4 - N + 3j 2 , 

C2 Vc 6 c 5 c 6 / c 4 c 5 c 2 



o 2 J4 J5 /C 2 C 3 C 2 \ ? 2 J 4 J 3 

3m — = — + 2— c 4 c 2 H (c 4 - c 5 ) 

C 4 Z \C5C6 Cg/ c 2 c 4 c 3 

j 5 2c 4 - 2c 3 - c 2 

+ ir v + 3j4 - 

2 c 5 I 
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(2) V0111 

3m 2 — = — (c 4 - c 5 ) + 2— (c 3 - c 4 ) - TV + 3j 3 , 
C3 c 3 c 5 

2m — = -— (c 4 - c 5 ) + -— (c 4 - c 3 ) H c 4 

c 4 3c 3 3c 5 c 2 

--7V + 2j 4 + j 2 , 

/, r - \ h ( . \ I r> «^*4 . .C 5 -C 3 C 2 + C 4 

2(N — 2j 2 -j4) = — (c 2 + c 4 ) + 2— c 2 + j 5 . 

c 5 c 4 c 5 c 6 

(3) V1011 

2m 2 - = ^ (c 3 - c 4 ) - iV + 2j 3 + j 5 , 

c 3 c 5 

(7V-2j 5 -i 3 ) = -(c 5 -c 4 ). 



c 3 



(4) V1010 



2m 2 3 J. = 3 -l CA — 2N + j2 + 2j 3 + 2j 4 + 2j 5 , 
c 3 c 2 

m 2 — = — (c 3 - c 4 ) - — c 4 + TV - 2j 4 - j 2 - j 5 , 

C5 C5 c 2 
, 2 J2 2 J5 C 2 j 2 . C 2 

4m —=m c 4 + (j 5 + 2j 4 ) — 

c 2 c 5 c 4 c 2 c 4 

-^-(c 3 -c 4 )-7V + 3j 2 , 
c 5 c 4 

4— (j 2 + 2j 4 ) = -m 2 ^ + 3 -* (4c 3 - 3c 4 - 4c 5 ) - (j 5 + 2j 4 ) ^ 
c 4 c 5 c 4 c 2 c 4 

+4 (j 2 + 2j 4 ) + (c 3 - c 4 ) 

c 4 c 5 c 4 

+9iV - 16j 5 - 8j 4 - 7j 2 , 

2j 4 + j 2 - 2j 5 = 2— (c 5 - c 4 ) - — c 4 . 

c 3 c 2 



(5) V0110 



m 2 — = — (c 3 - c 4 ) - — c 4 + TV - 2j 4 - j 2 - j 5 , 

c 2 c 5 c 2 

- 2 J5 2 h c 5 — C 3 ? 2 j 5 

4m — = m 1 (c 5 — c 3 ) H (c 3 - c 4 ) 

c 5 c 2 c 4 c 2 c 5 

+ (j2 + 2j 4 )^^-iV + 3j 5 , 
c 4 

, m2 /- ioM , 2 isC 3 2J2C3-C5 . J2, v ..C2 

4 — (j 5 + 2j 4 ) = 4m Ym V — (c 3 - c 5 ) - 4j 5 — 

c 4 c 5 c 4 c 2 c 4 c 2 c 4 
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+32 2— (4c 2 - c 3 + c 5 ) 

C4 C4 

+ ^("4^_4 C2 + 3 C 3_3 C4 ") 
c 5 V c 4 / 

+9iV - 16j 2 - 8j 4 - 7j 5 . 

(6) V1001 

Js_ (N-2j 3 -2j 5 ) (N-j 3 -j 5 -1) 
c 3 (N-2j 3 -2)c 6 

3m 2 ^ = 2^c 4 - f ^ + J ±) c 2 - N + 3j, + 3j 6 , 
3m 2 ^ = 2^ + ^c 2 -^c 4 -iV + 3j 2 . 



C 2 VC4 CqJ c 2 

(7) VOOIO 

m 2 33 = 3_3 ^ _ c ^ _ N + + 

c 3 c 3 

m 2 — = 2— (c 4 - c 5 ) + — (c 3 — c 4 ) + iV — 3j 5 , 
C5 c 3 c 5 

ra 2 — = 2— (c 5 - c 4 ) - — c 4 + 2j 5 - 2j 4 - j 2 , 

C2 c 3 c 2 

— (AT + 2 J4 - 2j 5 ) = 2m 2 ^ - (2 J4 + J 5 ) - + ^- (c 3 - c 4 ) 

C4 C 3 C4 C4 C5 C4 

_ m 2 J5 C2 + _ _ _ _ 

c 5 c 4 

(8) VOOOl 

j 3 _ (N-2j 3 -2j 5 )(N-j 3 -j 5 -l) 
c 3 (iV-2j 3 -2)c 6 

N - 2j 2 - j 4 = — c 2 - — - 2 - J 6 



C4 V 2 / cq 



C Analytical results 



For completeness, we present here the analytical results for all two-loop integrals shown 
in Fig.l 



Flllll(l, 1, 1, 1, 1, m) = -C(3) + \-^S 2 + 0(e), 
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F01111(l, 


1,1, 


1,1, m) 


= |c(2) In 3 - |C(3) + 3^=S 2 + 0(e), 

3 3 y/3 


FllllOfl, 


1,1, 


1,1, m) 


= --C(2) ln3 + -C(3) + 6^S 2 + 0(e), 


F00111(l, 


1,1, 


1,1, m) 


= 9-^5 2 + <9(e), 


F10101(l, 


1,1, 

I I 


1, 1, rn 

II / 


= -4C(3) + — -^=S 2 + i-C(2) + 0(e), 


F10110(l, 


1,1, 


1,1, m) 


= -C(3) + 9^5 2 + 0(e), 
V3 


F01101(l, 


1,1, 


1,1, m) 


= 6C(2) ln2-^C(3) + 0(e), 


F01100(l, 


1,1, 


1, 1, m) 


= ^^=S 2 + iir((2) + 0(e), 


FOOlOlfl, 


1,1, 


1,1, m) 


= -3C(3) + — ^=S 2 + ivrC(2) + (9(e), 


F10100(l, 


1,1, 


1,1, m) 


= -2C(3) + 9-^S 2 + ^ttC(2) + 0(e), 


F00110(l, 


1,1, 


1,1, m) 


= -iLi-2C(2)-2C(3) + 0( e ), 


F00100(l, 
F00001(l, 
FOOOOO(l, 


1,1, 
1,1, 
1,1, 


1,1, m) 
1,1, m) 
1,1, m) 


= -3C(3)-2i7rC(2) + C(e), 
= _3 C (3) +CTC (2) + 0(e), 
= -6C(3) + 0(e), 



where 



F{.ABJJ7C}(a, 6, i, j, fc, m) = m 2 K" 2 J d N hd N k 2 

P {a \k l ,Am)P (b) {k 2 ,Bm) 
P {i) (k 1 -p,lm)P U) (k 2 -p,J 
p( k \ki-k 2 ,lCm) 



9 9 ' 



and 

zr= r( l +£) , pW(fc,m)= , l2 1 

(47r)^(m2) e (P + m 2 )< 

and where the normalization factor l/(2n) N for each loop is assumed; A, B,X, J , K, — 0, 1 
and 

5 2 = ^=C7 2 (? ) = 0.2604341376 
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+- ? = In 3 + e| y - 6C(2) - -((3) - \2-j= - 63S 2 + 4£(2) In 3 + 9S 2 -^= 
A 2 In 3 + 4^ In 3- i^ln 2 3 - H_L C(2) _ Hl£fiM \ + ^ 
V0111(l,l,l.l.m) = ^ + ^ + f-|c(2)-^-^ 

+3^=113-9— li^|+0( £ 2 ), 

™*^-s + ;(!-5i) + T^^-« 

1 5 19 
V1110(l, 1, 1, 1, m) = — + - + - - 4C(2) 

+e{y - 20C(2) - 14C(3) + 24((2) In2| + 0(e 2 ), 

1 S 19 7T 

V1010(l,l,l,l^) = - + - + T + C (2)-3 7 = 

+ 2C(2) - C(3) - 15-^ - yS 2 + 9-^ ln3 j + 0(e 2 ), 
V0110(l, 1, 1, 1, m) = ^ + J (J - iTr) + y - 5C(2) - 3-^= - 2m 
+e| y - 10C(2) - C(3) - 15-^ - y S 2 + 9-^ In 3 + 2m [C(2) - 2] ) + 0(5 2 ), 



V 1 00 1(1 , 1 ,U, m )^ + ig-^) + f + 5 C(2) - 4 ^-f S2 
+-| In 3 + fjy + 8C(2) + §C{3) - 12-?= - 63S 2 + ^S 2 k3 

V0011(l, 1, 1, l,m) = ^ + A + y - 2C(2) + 4f " 6C(2) - 4C(3)| + 0(e 2 ), 
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V0010(l, 1, 1, 1, m) = — + - (- - iir) + — - 7C(2) - 
+e! { j - 17C(2) - 11C(3) + iTT [2C(2) - 7]| + 0(e 2 ), 

V1000(l, 1, 1, 1, m) = ± + A + ^ + 2C(2) + ejy + 10C(2) + 4C(3)| + G(e 2 ), 

1 5 19 f 65 1 

V0001(l, 1, 1, 1, m) = — + - + - + C(2) - in + e| - - 3((2) - C(3) - 7m j 

1 1/5 \ 19 
V0000(1, 1, 1, 1, m) = — + - - zttJ + — - 7C(2) - 5m 

+e[y - 35C(2) - 5C(3) + «r [6C(2) - 19] | + 0(£ 2 ), 



V{IJICC}(i,j,a,b,m)=K- 2 J d N k l( i N k 2 P^(k 2 - p,lm) 

P {3 \k x -k 2 ,Jm) 
P {a \k u Am) p( b \k 2 ,Bm) 



and 



Ls, ( y ) = -2.14476721256949 



where we use the following definition [27]: 



Ls 3 (x) 



In 2 



6 

2 sin - 

2 



d0. 



To our knowledge LS3 yf-J has appeared for the first time in the calculation of the £-part 
of two-loop tadpole integrals in Ref. [ 



17 59 



65 



Jlll(l, 1, 1, m) = -m 2 — + — + — + e { — + 8C(2) 



2e 2 4e 8 
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1117 



-e- 1 — - 52C(2) + 48C(2) In 2 - 28C(3) j + 0(e »)j , 
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7T 



+14 7! c(2) + 27 Vs 



Lss{ 2 f) 3 



J011(l,l,l,m) 



-^C(4)| + 0(e 3 ), 
11 35 3 7T 



5 39 7T 

-m" I — H 1 1 = + £< 1 = 

U 2 4e 8 2^3 \ 16 4^ 



9 vr 

"273 

117 7T 



1 



} f 1033 _ 1053 n 345 vr , 243 
|~32 



-i — =ln3 - ^^ln 2 3 - 21-^=C(2) - 



4 ^3 



27 7T 

T7t 



J001(l,l,l,l,m) 



8 y/3 ' 4 
1 5 11 f 55 



-m 



2e 2 4e 



16 



<^ 3 )J, 
-5C(2)-4C(3) 



f 949 1 1 ~1 

-e 2 ^ - yC(2) - 10C(3) - 32C(4)| + 0(e 3 ) 



J\ 13 m [115 13. 

JOOO l,l,l,l,m = m 2 — + + e<^ ( 2 zvr 

Ue 8 2 I 16 2 SV ; 4 



865 91 



32 4 LC(2) -2 C(3)+ - 



3C(2) 



115 



0(e 3 



where 



3{lJC}{i, j, I, m) = R- 2 J d N kid N k 2 P (i) (k 2 , 2m) 

P®(ki - k 2 , Jm)P m {h - p, Cm) \ p2 __ 



and for the £ 2 -part of the J011(l, 1, 2, m) integral we take into account the "guessed" 
expansion from Ref.[26]. The expansion of Jlll(l,l,l,m) up to e A terms is given in Ref.[29]. 
A very elegant e-expansion up to arbitrary order for the one-loop propagator type integral 
can be found in Ref.[30]. 
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